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The Hanbury Brown-Twiss (HBT) effect in two-particle correlations is a fundamental wave phe- 
nomenon that occurs at the sensitive elements of detectors; it is one of the few processes in elementary 
particle detection that depends on the wave mechanics of the produced particles. We analyze here, 
within a quantum mechanical framework for computing correlations among high-energy particles, 
how particle detectors produce the HBT effect. We focus on the role played by the wave functions of 
particles created in collisions and the sensitivity of the HBT effect to the arrival times of pairs at the 
, detectors, and show that the two detector elements give an enhanced signal when the single-particle 

wave functions of the detected particles overlap at both elements within the characteristic atomic 
transition time of the elements. The measured pair correlation function is reduced when the delay 
in arrival times between pairs at the detectors is of order of or larger than the transition time. 

PACS: 25.75.-q, 25.75.Gz, 03.75.-b 
, Keywords: Particle interferometry, GGLP effect, HBT effect, Relativistic heavy ion collisions 
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I. INTRODUCTION 
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The Hanbury Brown-Twiss (HBT) effect in nuclear and particle physics PH| is the enhancement at small relative 
momentum in the probability for observing pairs of identical bosons drawn from the same collision event, compared 
to that for pairs chosen from different events. Underlying the correlations between the fluctuations at two nearby 
but separate detectors is the symmetry of the wave function of a system of identical bosons under interchange of any 
two particles. HBT particle interferometry is in a special class of high energy physics experiments that, like kaon 
regeneration and neutrino oscillation searches, study an effect basic to quantum mechanics: amplitude interference. 
However, while, e.g., kaon regeneration experiments examine how weak- interactions affect the internal degrees of 
freedom of kaons, HBT interferometry probes the many-body spatial wave function of the final state particles. The 
effect, as we emphasize here, is fundamentally a wave phenomenon manifested at the detectors. Indeed, as is well 
understood in terms of classical wave mechanics, HBT intensity-intensity correlations may be observed with any type 
of wave satisfying the superposition principle (e.g., sound waves) with at least two independent incoherent sources 
and two detectors. 

In this paper we analyze from a quantum mechanical point of view the physics of correlations between identical 
particles, and how the detectors in high-energy physics experiments probe the many-particle wave functions in space 
and time. While the spatial requirements on the detector separations for measuring an HBT effect are well understood, 
the temporal requirements are less well studied. We may ask, what is the maximum separation in time of arrival of 
two identical particles at the detectors that will still yield an HBT effect between the two particles? It is commonly 
assumed in collision experiments that interference happens only between particles produced in the same event, not 
between particles from different events. One could imagine that this occurs because the correlations are built in at 
the time of the event; however since one detects interferometry between photons from the opposite edges of a star, 
where there can clearly be no correlation between the emission processes, this is certainly not the case. Underlying 
the assumption of lack of interferometry between particles from different events is the existence of a timescale for 
detection of interferometry, which is exceeded by successive events. As we show here, the crucial timescale for HBT 
correlations is that inherent in the detection process; the maximum value of the measured HBT pair correlation 
function is reduced when pairs of waves arrive at the detector atoms with delay times comparable to or larger than 
the characteristic transition times of the atoms in the detector. 

Considering the wave mechanics of the emitted particles in space and time also enables us to understand more 
fully questions such as: How does the exchange symmetry of the many-particle wave function lead to detection of 
momentum correlations? How does the wave packet nature of the single-particle states affect measurements of the 
HBT effect? For example, how docs the spreading of wave packets affect the measurement of correlations? 
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Part of the motivation to examine the role of the interplay of the time scales of the pion field and the detection 
process in particle interferometry experiments is the clear illustration in the original studies of photon correlations 
of how detection times affect the size of the maximum value of the measured two-photon correlation function • 
As Purcell || pointed out, for chaotic stationary sources the maximum height of this function above unity is ap- 
proximately the ratio of the coherence time of the light to the observational binning time. Hanbury Brown and 
Twiss's measurements on the star Sirius JlO[ ] produced a maximum correlation signal ~ 10~ 6 , in contrast to data 
on like-pion pair correlations gathered to date, with final-state interaction corrections included, which indicate that 
the two-pion correlation function, C2 ( 1 9* I ) 5 rises approximately (but not unambiguously) to two in the lowest bin in 
relative momentum. 

Four time scales are basic in a measurement of the HBT effect. The first is the lifetime of the source, r s , over which 
it emits particles; in a high-energy collision the production region exists for ~ 1-10 fm/c [~ 0.3 — 3 x 10 -23 sec] in its 
own rest frame. Second is the coherence time of the source, r c , the characteristic time for each elementary radiator 
to emit a wave. Generally, the formation time for a pion or other particle depends on its energy; the characteristic 
lifetime of elementary pion sources r c is < 1 fm/c, and t c <t s . The third scale is the atomic transition time, r a tomic, of 
the energy-absorbing material of the detector, the time over which a mobile electric charge is created in the detector; 
on this time scale the detector atoms "do quantum mechanics" on the incoming particle waves, i.e., the atoms are 
sensitive to the amplitudes and phases of the waves. For ionization of a gas atom, T atom ic ~ h/10eV ~ 10 _16 sec. 
The final scale is the exposure time, r cxp , the interval between observations of the state of the detector atoms. The 
time of an accurate momentum measurement on a relativistic charged particle in a magnetic spectrometer of length 
^spec is ~ ^spec/c, which for a typical length scale of 10 m gives a characteristic time to measure the momentum 
7"cxp ~ 3 x 10 -8 sec. A more familiar but less relevant measurement time scale is the detector resolution time, r rcs , 
literally, the minimum time to detect an electronic signal, which includes the time it takes to collect and amplify 
the initial electric charge. A typical resolution time for a wire chamber, essentially the rise time of the voltage pulse 
produced by an electron avalanche on an anode wire, is on a nanosecond scale, r rcs ~ 10~ 9 sec. As we show in Sect. Ill, 
^atomic, rather than r rcs , is the important time scale over which HBT correlations are detected. 

This paper is organized as follows: In Section II, we develop a framework for understanding the temporal structure 
in correlations between quantum mechanical particles created in high-energy experiments. To be specific, we describe 
pions, but our results hold for other particles as well. In Section III, we analyze the features of relativistic wave packets 
important to HBT, and compute the detection probability for pairs directly-produced in a heavy-ion collision and 
then for those from direct production plus resonance decay. In Section IV, we summarize our results and conclusions. 
In the Appendices, we describe the details of the quantum mechanics of particle detection. 



II. DESCRIPTION OF THE PROBLEM 



One and two pion measurements on a multiparticle system are described by the single-pion and two-pion density 
matrices for particles of given charge: (<fi^ (xi)4>(x2)} and {<j>' {x\)(f^ {x^<j>{x^)(f){x^)) , where 4>(x) is the part of the 
(Heisenberg representation) pion field operator that destroys particles of the given charge, and the brackets indicate 
an ensemble average over the states of the colliding nuclei. The pion-pair correlation function, C 2 {q), which depends 
directly on these two functions, is measured as the ratio of the pion pair distribution to the separate single distributions: 

C(n\= {d 6 n 2 /dp 3 dp' 3 } 
2[q) {d 3 m/dp 3 d 3 ni /dp' 3 Y y ' 

where q = (p — p')/2. The braces in the numerator denote an average over an ensemble of pairs drawn from the 
same event and in the denominator they denote an average over an ensemble of pairs drawn from different events. 
The single-pion momentum distribution d 3 n\/dp 3 is given in terms of the plane- wave momentum state creation and 
annihilation operators at and a ff by 



^ = (4a-> = 2e p jd 3 ri d 3 r 2 e -*H*x-* 2 ) (tf fafrfa)) , (2) 



with e p = (p 2 + m 2 ) 1 / 2 and p ■ x = Lot — p ■ r, where w is the particle energy and p is the momentum. Note that the 
explicit dependence on t\ and t 2 in the phase factor of Eq. (^|) is canceled by the time-dependence of the single-pion 
density matrix, as later shown in Eq. ( |l0| ) . The momentum distribution of pairs is given similarly by 
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^jSs = (44' V a ?) = A£ p £ p' J dVl ' ' e-^-(—^e-^'-(— 3 ) (^(a!i)^(x a )0(x3)^(x4)) ■ (3) 



A detailed description of the pion wave functions that emerge from the production region of a heavy-ion collision 
requires knowledge of the evolution and geometry of the source and how field excitations of quarks, gluons, nucleons, 
etc. form the currents that radiate pions. However, these details are unimportant for our present purpose of under- 
standing how we detect HBT. We retain only the essential features of the pion production process, in particular its 
characteristic length and time scales. In high-energy collisions, pions are produced in bremsstrahlung-likc processes, 
in states similar to momentum wave packet states. Final-state interactions, particularly Coulomb and strong interac- 
tions, severely influence particle state evolution; also, collisions between pions and air molecules, detector materials, 
and other target nuclei can significantly alter C 2 ] in particular, multiple scatterings of particles after they leave the 
source redistribute pion tracks thereby reducing the size of the measured correlation in each bin in relative momentum 
and decreasing the effective size of the source measured in HBT experiments (l^Jl^] . Thus in order to focus on the 
wave mechanics of detecting HBT, we assume here that after their last strong interaction the wave packets propagate 
in vacuum. 

The production of charged pions and their propagation to the detectors is described by the Klein-Gordon equation 
for the charged pseudoscalar field, which relates the pion field, 0, to the source, J(x), of the pion field at the last 
strong interaction: 

(d 2 /dt 2 - V 2 + to 2 ) <P{x) = -J(x), (4) 
where is the part of the current operator that emits pions of the given charge. Thus 

<P(x) = { dx' D Iet (x - x')J(x') , (5) 



where the free-field retarded Green's function, D let (x — x'), satisfies 

(d 2 /dt 2 - V 2 + to 2 ) D Tet (x - x 1 ) = -5 {i \x - x 1 ), (6) 
and vanishes for t <t'\ the integrations are over all space and time. The Green's function has the representation 

d 4 k e -*k-(x~x') 

(2^)4 ( w + ie) 2 - e\ 



D Ict (x -x')= / t^3j 7 ■ ■ .^ 2 _ 2 , (7) 



where e is a positive infinitesimal and u> — k . Integrating over ui for t > t', we find that the created pion field is given 
by 

^ =/^SrF L dt i d3r ' jMe " ,HW ' (s) 

with k on-shell. Since one measures at t much larger than the source lifetime we extend the upper limit in the time 
integral in Eq. (^J) to +00 and write 

/d 3 k 

where J(k) is the Fourier transform in space and time of the pion source current operator. Equations ^ and (||) 
imply that the Lorentz-invariant one pion momentum distribution is given by 

2 £p ^-= [dx 1 dx2e- ip - (xi - X2 \j\x 1 )J(x 2 )) (10) 
dp J 

and the two pion momentum distribution by 

4£pgp> = / dxidx 2 dX3d Xi e -iP-(*i-*4)-iP '-(xa-ss) (j\xi)j\x 2 )J{x 3 )J{x 4 ))- (11) 

dp dp' J 



Eqs. (JlOj) and ( |ll| ) show how the HBT correlation function directly probes the correlation functions of the source, a 
point of view initially introduced in Ref. B . 
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The HBT effect arises only when single particles are produced in mixed quantum states, that is, in an ensemble or 
statistical mixture of single-particle states. (See Ref. jl3| for further discussion of the multiparticle states in HBT.) 
When particles are produced completely independently, e.g., by a thermal source, the source currents factorize as 

(JHx 1 )JHx 2 )J(x 3 )J(x 4 )) = (JHx 1 )J(x 4 ))(J\x 2 )J(x 3 )) + (j\x 1 )J(x 3 ))(J^x 2 )J(x 4 )). (12) 

The pion correlation functions similarly factorize, and the momentum distribution of pion pairs |D| is 

d 6 n 2 d 3 n 1 d 3 n 1 1 



dp 6 dp 1 ' 6 dp 6 dp' 3 ' 4ep-e ? , 



\(j\p)J(p'))\-i (13) 



HBT interferometry seeks to measure the second term in Eq. (|T^). The key physical mechanism that leads to the 
factorization of currents is the loss of phase correlations among the elementary sources, which is expected to occur 
in heavy-ion collisions through the considerable rescattering of pions in the production region. [] The HBT effect is 
maximum, for incoherent emission. 

By contrast, when particles are produced completely coherently, e.g., as in an atom laser beam extracted from a 
Bose-Einstein condensate or by an ideal chiral condensate |l7j , the source currents factorize as 

{J\x 1 )J\x 2 )J{x 3 )J{x 4 )) = (J^x 1 ))(J^x 2 ))(J(x 3 ))(J{x 4 )) . (14) 

In this case, single particles are produced in a pure quantum state; the momentum distribution of pion pairs is the 
product of two single-pion momentum distributions, and the HBT effect would be absent. 

III. MEASURING CORRELATION FUNCTIONS 

We turn now to the physics of momentum measurements. How ionization chambers (and many other detectors used 
in high-energy physics) track well-localized particles is understood from semi-classical ideas; however, the response of 
such detectors to a particle whose wave function is not well-localized in space and time necessarily requires a quantum 
mechanical description (a problem of longstanding interest, e.g., from Ref. Jl8| to (l9).) HBT correlations are the 
direct result of the interactions of the many-pion wave function with the electrons of the energy absorbing material 
of the detector. Computing the probability of measuring the momentum of a fast charged particle is a problem in 
multiple scattering, where we must compute the probability of measuring a particle track through a spectrometer 
system. The first interaction selects out the direction of the momentum, and subsequent interactions with more atoms 
and a magnetic field select out the magnitude of the momentum. We define the single-pion momentum measurement 
probability P^(k ) as the probability for a pion to ionize a detector gas atom at some location a, say the first atom along 
a track, and undergo a transition to a plane- wave momentum state k. Similarly, we define the pion-pair momentum 
measurement probability P£™ {k ,k' ) as the probability for each pion to ionize one atom and emerge in a plane wave 
state, one atom at a with a pion in the final state k and one atom at b with a pion in the final state k' . 

As derived in Appendix A, the crucial function describing the response |2(| of the detector in measuring a particle 
of momentum k - a "momentum measurement" as defined above - is the spectrometer function, 

S^(xix 2 ) = e ^{xi) (j a (x-L)j a (x 2 ))il)^(x 2 ), (15) 

where j a (x) is the effective electromagnetic current operator for the atomic electrons, Eq. (pl|), and ip^ is the final pion 
plane- wave momentum state wave function. As shown by detailed calculation in Appendices A and B, the momentum 
measurement probabilities are given by the pion correlation function as filtered by the "spectrometer function," i.e., 
the overlap in space and time of the single and two-pion correlation functions (for given final pion states) with the 
correlation function of the effective electron currents in the detector atoms, 

P^(k) = \dx x dx 2 SUx lX2 ) (tffaMxi)), (16) 



1 In the absence of rescattering the statistics of the constituents of the currents can destroy this factorization; a simple example 
is the correlations among pions radiated by a weakly- interacting gas of nucleons |ll|,|l4| ; see also [|L3| . 



4 



and 

P£(k,k')= I dx 1 dx2dx 3 dx i SUx 2 x 3 ) (^{ Xl )^(x 2 )cj)(x 3 )^(x 4 )) Sl(x lXi ). (17) 



When the source term factorizes as in Eq. (|l2|), we obtain the pair momentum detection probability as a sum of 
two terms, 

P£(k ,k') = D*{k ,k')+E%(k,k'), (18) 

where the direct term is 

D*(k,k')= [dx 1 dx 2 dx 3 dxiS%(x2X3)Sl(x 1 Xi) (^(x 2 )^(x 3 ))(^(x 1 )<j)(x i )), (19) 
and the exchange term is given by 

E^(k,k') = f dx x dx 2 dx 3 dx A SI(x 2 x 3 )Sl(x lX4 ) (^(x 1 )cf ) (x 3 ))(^(x 2 )^(x i )) . (20) 



The direct term is the product of the single- pion detection probabilities, Eq. (16), 

Dt(k,k') = p:(k)p ] :(k'), (21) 

and is equivalent to the probability of observing two like-charge pions from two different collision events, one in state 
A; at a and one in state k' at b . The exchange term is the HBT effect. 

Knowledge of the mixture of pion states produced by the source is critical to understanding HBT interferometry. 
For pions produced in single particle states <fi{x) the single-pion correlation function is 

(^(xt)^)} =J2 F i ¥>*(*!)¥>,• (* 2 ), (22) 

i 

where Fi specifies how the distribution of pion states depends on the evolving geometry of the particle production 
region. [Generally, the decomposition of the single pion correlation function as a sum of single particle states defines 
the single particle states.] The corresponding pion source current-current correlation function is 

(ji( Xl )J(x 2 ))=J2 F i JfMJifa), (23) 

i 

where J~i(x) is the transition matrix element of the pion source operator J(x) . To illustrate the structure concretely, 
we approximate the sources of the pions produced in a high-energy collision as Gaussians in space and time 

J.( X ) = f ^ P -(*- Xo ) e -(t-to) 2 /2r? e -(f-f Q f/2B^ (24) 

(2ir) t c R c 

where R c and r c are the characteristic length and time scales for pion formation, and N is a normalization constant; 
here the subscript i stands for the central momentum p and the space-time origin of the wave Xq. It has been shown 
elsewhere [p p|j2l| that the effective relative momentum scale of the pair correlation function involves not only the 
dimensions and lifetime of the particle production region and the mixture of states produced but also the individual 
particle formation length and time scales. These sources give rise to pions in Gaussian wave packets, which in the far 
field have the form 

1 r°° 

P<(*)»-o-5I / de q J i (e q ,qf)e-^ t+i " r , (25) 

where Ji (e q ,qf) is the four-dimensional Fourier transform of Ji [x] . Equation (^j|) follows from Eq. (^) by selecting 
the outgoing wave after integrating over momentum directions. 

For p = pz, the Gaussian source gives rise to a wave packet which spreads out approximately as 
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where j p = Sp/m^ and v p = \p\/e p . For t ^> t c , the ratio of the longitudinal to the transverse spread of the wave 
function is ~ l/7p • Far from the source, relativistic pion momentum wave packets are pancake-shaped and move in a 
direction perpendicular to the face of the pancake. Consider, e.g., a 1 GeV pion and assume that r c ~ R c ~ 1 fm. For 
a source-detector distance L~lm the size of the wave packet envelope at the detectors in the transverse direction is 
~ 20 cm and is ~ 0.2 cm thick in the longitudinal direction. 

Substituting Eq. (||) into Eqs. (|l|) and @, we find the direct and exchange terms, 

B* (k k' ) = £ FiFj Wl(ii)* W| (jj) (28) 

ij 

Et{kk' ) = J2 F i F j Wl (ij) . (29) 

ij 

where we define the wave function-detector overlap function 

Wl(ij) = jdxdx' i p*(x)SI(xx') 9j (x'). (30) 

Each possible pair of pion waves (ij) contributes to the direct term, with one detected at a and the other at b. 
However, only pairs of waves that arrive together at each detector atom, within the time interval set by the detector 
functions S contribute to the exchange term. Interferometry occurs when the wave functions of each of the particles 
overlap in each of the detectors at the same times. Equation (|29|), with (psj), shows that there is no restriction on the 
time interval between the arrivals at the two atoms. Correlations are measured over the time scale imposed by the 
response of the atoms, Eq. (|3^), the time scale over which the atoms are sensitive to the phase and amplitude of the 
incoming particle waves. 

We illustrate the physics of detecting correlations with a simplified model of the effective atomic current-current 
correlation function in Eq. (|l5|), which takes into account the main features of the dynamic response of atoms to 
ionization by relativistic charged particles , namely, we assume that the space and time dependence of the effective 
atomic current-current correlation function factorizes as 

(jl(x 1 )j a (x 2 )) = f a (r 1 f 2 )g(t 1 t 2 ). (31) 

Since the statistical distribution of atomic electron states in a gas under normal conditions is essentially time- 
independent, then g(t\ t 2 ) = g(t\ — t 2 ) . The Fourier transform of g(t) is the atomic energy-absorption spectrum, and 
is determined from the distribution for energy-loss per ionization for a relativistic charged particle passing through a 
gas. In a monatomic gas, this distribution is approximately Gaussian with a long tail extending up to the kinematic 
limit for energy transfer, 2m e (/?7) 2 . The majority of ionizing collisions occur within the Gaussian part of the distri- 
bution and it is particularly the ionization events in which the free electron carries away a minimal amount of kinetic 
energy that are important for tracking. Thus, we can neglect the high-energy tail and select a (normalized) Gaussian 
energy spectrum for g: 

9(t)= l^ffHe--*, (32) 

where g(u>) — -\/2n C -1 e _ ^ _CJ0 ^ / 2c = . The average energy loss, ujq , for creating electron-ion pairs in, e.g., noble gases, 
is 10 — 50 eV H^]. The energy bandwidth of the Gaussian part of the distribution, £, is ~ 10 eV, an energy on the 
scale of the average ionization potential per electron in a Thomas-Fermi model of an atom. This scale determines 
the characteristic time for ionization, T at0 mic = 1/C ~ 10~ 16 sec. The spatial scale of the atomic correlation function 
is determined by the sizes of the atomic electron wave functions, characteristically ~ 1 A; in terms of momentum, 
if the kinetic energy picked up by an electron is a few tens of eV then the electron momentum, \q \ , is on a keV/c 
scale, corresponding to a distance l/|<f| ~ 1 A, a result roughly consistent with the size of impact parameters required 
in a classical Weizsacker- Williams picture of a collision. For convenience, we model f a (r\ r 2 ) as the product of two 
Gaussian functions centered on the atomic nucleus: 

/ a (rl f 2 ) = f a e -{r 2 -af/2Rl (33) 

where i? a is of order angstroms and / a is a constant. 
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The delays between the emissions of pions directly produced in a high-energy nucleus-nucleus collision are no larger 
than the lifetime of the source. However those produced in resonances can have a longer spread in emission times. As 
we now show, Hanbury Brown-Twiss correlations are insensitive to the delays in arrival times at the detector atoms 
generated by emission delays much smaller than the atomic response time, T a t om ic- With the use of the far- field form 
of the wave functions for directly produced pions, Eq. with Eq. the space and time integrations in Eq. ( |30| ) 
are trivial. Integrating over the energy components of both pion wave packets we find the intermediate form, 

WI(ij) = c a J^~g(u)e- i ^+^o-t' ) 

x e -[k- q afRl e W-<f a -rS)j*( ek + UiqSL ) J ff ,(e k +u;,qa), (34) 

with c a = ne 4 f a R®/ 4a 2 EkV ; \q \ is determined by the condition e q — Sk + u), and t — t' is the emission delay time 
between the pion waves. Since as one sees from the singles distribution |23]] , the transition matrix elements vary over 
an MeV scale or more, and g(uS) restricts w to a neighborhood of size £ about ujo , both on an eV scale, then for 
pions of energy at least one GeV, it is a very good approximation to replace q with k and neglect to everywhere in the 
second line of Eq. (^) . Consequently, the first term on the second line of Eq. ([34]) requires that fc • a w 1 (to within 
one part in at least 10 10 ), so that k = ka . The integral over u> is then simply the Fourier transform of g : 

g(t - t' ) = e -^o(t -t' ) e -(t -t' fC 2 /2 ^ ^ 



« c a g(t - t' ) J2 x (e k , k ) J , x , (s k , k ) . (36) 



and the overlap function is 



The direct term is independent of g, since t = t' and p = p' in Eq. (p6|). In fact, Eqs. (|3J]) and (|3(]|) show that when 
the emission delay time between pions is much less than 1/C = ^atomic we may also neglect the time dependence of the 
detection process in the exchange term; to a very good approximation we may then write \g(t — to') I ~ l<?(0)| = 1 • 



Substituting Eq. (|36|) into Eqs. ([28|) and (29) and summing over the mixture of wave packet states we find that the 
pair momentum detection probability is proportional to the Fourier transform of Eq. ( p^ ) , 

P£(k,k') = c a c b [{jHk)J(k)){jHk')J(k')) +|<J t (fc)J(fc')>| 2 ] , (37) 

the naive result obtained by neglecting the time dependence of the detection process. However, in situations leading to 
much longer time delays, in particular, pion emission from long-lived resonances compared with direct pion production, 
one cannot necessarily neglect the time dependence of the detection process.^ 

Consider interferometry between a 7r~ produced directly in a heavy-ion collision and from the decay of, say a 
lambda, A — > tt~ + p, produced in the same reaction. Because the A moves more slowly than a directly produced 
pion of the same rapidity as the one emitted in the decay, the pion from decay will lag the directly produced one 
by a time 5 at the detector atoms. To estimate this arrival time delay, we note that a 7r~ emitted in the forward 
direction has rapidity yo ~ 0.67 in the A rest frame, and that a A of rapidity y travels on average a distance t\ sinhy 
before decaying, where ta is the A lifetime. Thus, S = T\/(coshy + sinhy/ tanhyo); which for a A of typical rapidity 
3 is ~ 0.037ta = 9.7 x 10~ 12 sec, much longer than the atomic time scale. Pions emitted in other than the forward 
direction will have an even greater time lag. As we shall see, detector atoms are sensitive to such delays when we 
measure pair correlations. 

The reduction in the probability for detecting pairs due to time delays between direct and resonance decay pions 
can be readily estimated using a simplified scalar-field model to compute the overlap, Eq. (j^Tj), with a detector atom of 
the wave functions from a direct-production pion and a pion from, say, lambda decay. The wave mechanical features 
of the exact problem do not depend on the details of the model. Consider the interaction Tlj(x) — a n' (x)p(x)A(x), 
where a is the coupling constant. The transition matrix element for the decay-product pion is 

jf(x) = ar p (x)<P A (x), (38) 



2 Effects of time delays in propagation of wave packets can also enter in the observation of neutrino oscillations, as noted by 
Kim [M, and references therein. 
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where ip A is the wave function of the A baryon of width T ~ 10~ 6 eV, and ip is a plane wave state for the proton with 
momentum p. As before, we write the wave functions of the lambda and both pions in the far-field form, Eq. J25|). 
The Fourier transform of Eq. ( |38| ) is 

7 A/_ ff) * J A (e + el,g+p) 

^TT\ t i'i)— . I / . p A \ I AT /O ' ' 

2 4 +P J^y {£ + £ p " £ ?+p ] + lT/2 



where Jk, the source function for lambda production, is similar in structure to Eq. (24), e~ = (q 2 + rrijj 1 ^ 2 , and 



eZ = [p 2 + Trip) 1 / 2 . We subsume the origin and momentum labels of the wave packets into the subscript on the 
transition matrix element, J7a- Thus, the wave function overlap at a detector atom located at a with pion final state 
k is 

fe( P) a K*j*& J toiu + el + Z-efo+W m 

The overlap function Wj? (it Ap) for an atom at b is exactly the same as above, but with k' in place of k. 
The contribution to the exchange term, given by 

P-( h A) (*,*')= E F - Fa W|(7r Ap )* H| (tt Ap), (41) 

decreases as the width T of the lambda decreases, i.e., as the lifetime of the lambda increases. We can see this 
explicitly by selecting a convenient form for the energy absorption spectrum of the atom: 

^ = (.-.o) 2 + (C/2) 2 ' (42) 

where ( ~ 10 eV (cf. Eq. fl3^)). The effect is simplest to see for k — k! (we report the calculation for k ^ k' in a later 
paper); then the sum in Eq. ( |4l| ) is 



~?* fA ^N 3 V ' (Hi/ (-o + £ l + 4-4 + P 2 + (C + r)V4- (43) 

Since the factor containing T varies over a much smaller energy scale than that of the rest of the integrand we can 
replace the last term in ( f43| ) with 27T(5(eI + eZ — e^ + ^)/(C + T) , where we have also neglected luq. It is convenient to 
assume that the A is produced in a spherically symmetric state in the lab, so that the Fourier transform of J A only 
depends on energy. The width of the decaying particle to lowest-order in our model interaction is T — a 2 po /87rm A , 
where po = ((mj/2mA) 2 — m 2 ,) 1 / 2 , and m 2 = m A + m 2 — m 2 . Thus, integrating over proton momentum we find 

P^f A) &,k) = y^K^A) $,k) (44) 



where 



2 P A 



V^f A) (k,k) = c a c b ^- £ F W F A \J: ( £ | , k ) J A (^ , q ) | 2 In | g ' ^ | , (45) 

^POK ff A tq - Ptt<70 

with q = Po m A/ m TT, P-kI-k — |fc|/ m 7rj an d replacing J A with an average matrix element. Equation ( [i"5"|) is the 
probability obtained when the atomic response time is neglected, that is, when one takes g(ui) = 2tt5(uj — u> a ). 
Equation (^) shows that pions from the decay of long-lived resonances can lead to arrival time delays at the detector 
atoms, relative to direct-production pions, that do indeed give a reduced contribution to the HBT signal. Thus, when 
r <C C the 7rA exchange term is reduced by a factor r/£ (for lambda-decay pions, ~ 10~ 7 ) from the naive result 
and for short-lived resonances with r>( we recover the case where we neglect the detector response time. Among 
the common weak decays that produce pions the detection time effect is strongest for interferometry with charged 
7r's from mesons with T < 10~ 7 eV; for AT ± , K^, T/( < 1CP 8 . We expect a weaker suppression for pairs from the 
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shorter-lived particles where T < 10 _5 eV; for Kg, A, and E , T/( < 1CP 6 . For rarer decays like S , f2~, and heavier 
flavors with T < 10 _3 eV the detection time suppression factors are less than 10~ 4 . 

The full exchange term, of which Eq. (^J) is a part, includes all possible pion pairs which involve contributions 
of the type: (J*,J V ,), (J*,J A ,), {Jt,J^), and (J£,J A ,). The {J*,J^) term is computed in Eq.©. The terms 
> ^A' ) an d G^A ' ^ir' ) are identical and contribute a factor 2 times Eq. (E3) to the full exchange term. The (J^ , J A , ) 
term may be computed in a way similar to that of Eq. (Q), giving an even smaller contribution to the HBT effect. 
The essential result is that the response time of the detector atoms is sensitive to phase differences between the waves 
due to arrival time delays at the detectors; a result independent of the origin of the waves. One can see this, quite 
simply, by introducing an explicit time delay in the wave funtions at the detectors in Eq.(|30|). 

IV. SUMMARY AND CONCLUSIONS 

We have developed, within an elementary quantum mechanical framework for computing correlation measurements 
in high-energy experiments, a general description of how detectors probe many-particle wave functions. The HBT 
effect is the consequence of wave mechanics performed by particle detectors and depends only on the wave functions of 
the particles at the sensitive elements of the detectors; it is not caused by stimulated emission or any other mechanism 
at the source and it does not depend on the history of the particles, e.g., the particles do not have to have a common 
origin. 

We have studied how momentum correlations between pairs of particles arc detected via the HBT effect. We have 
shown that: 1) The like-pair correlation function is able to reveal momentum correlations because the single-particle 
wave functions of the detected particles overlap at at least two sensitive elements of a detector, within the characteristic 
atomic transition times of those elements. 2) There is no restriction on the time interval between the transitions of the 
two detector atoms. 3) The size of the measured pair correlation function is reduced when the delay in arrival times 
between pairs at the detectors is of order of or larger than the transition time; e.g., delays from particles produced in 
very long-lived resonance decays. 
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APPENDIX A: MOMENTUM DETECTION PROBABILITY FOR SINGLE CHARGED PIONS 

We compute here the probability for an energetic charged pion to ionize a gas atom and emerge in a final plane- wave 
state with momentum k. We denote the initial state by \T) = \I, i a ,0), where / is the inital pion state, and \i a ) is 
the initial state of the atomic electrons; and the final state by \T) = \k, f a ) where \ f a ) contains the same number 
of electrons as |i a ), but with one electron in a continuum state. For simplicity, we work with single-pion states; the 
generalization to many-pion states is straightforward. The interaction of a charged pion with a detector-gas atom 
located at a space point a is 

H I =H 7T I +H^ (46) 

where H] — ej%A^, is the electromagnetic field, j% is the charged pion electromagnetic current, Ttj = ej^A^, 
with the electron current of the detector atom and — e the electron charge; we work in the interaction picture. For 

7r ± , the electromagnetic current is j%(x) — ±<^ (x)id^ (f>(x) . 

The lowest-order contribution to the amplitude for ionization in a collision between an energetic charged pion and 
an atom comes from the second order terms in the matrix element of the time evolution operator, U(t,t'), 

Al{T^T) = - f dx 1 dx 2 (J r \T[n a I (xi)n}(x 2 )}\l) = -ie 2 j dx 1 dx 2 {U^{x l )\QD^{x l x 2 ){k\fM\I). (47) 

Here T denotes time-ordering; the time integrations are from t to t' ', where the exposure time r cxp is t' — t, and the 
space integrals are over all space. In the latter expression we use the free-field forms for the operators, and introduce 
the photon propagator, D lxv {x\ x 2 ) = — i(0\T [A ll {xi)A u {x 2 )\ |0), where |0) is the vacuum. 
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The transition matrix element of the pion current becomes (k ^ {x 2 )id 2 (j>(x 2 )\I) = ^(x^id^ (0|<^(o; 2 )|/), where 
tp~(x) = (1/ ' \/2£kV)e~ %k ' x . Since the energy spectrum of the atomic states important for tracking via ionization only 
has Fourier components a few tens of eV or less, the spectrum of the atomic function J dxi(f a \jg(xi)\i a )D flv (xi, X2) 
is similarly constrained. Integrating in x 2 by parts in Eq. (47) and neglecting the derivatives of this function, we can 

replace id^ by 2k v . 

Squaring the transition amplitude, summing over all final electron states, averaging over initial pion and electron 
states, we find the ionization probability 

!»{%) = 4e 4 y dx!dx 2 dx3dx 4 £; CT (£3x4)0^3)^^ (48) 
where the single-pion density matrix is 

(0t(x 4 )0(x 2 )) = Y i Pl(lW{x*)<K*2)\I), (49) 

with pi the probability that the state |7) is produced by the source; the electric current-electric current correlation 
function for the atomic electrons is 

(£(s3)#(a0) ^E^E^IJa^)!^)^^^!)^), (50) 

/ 

where pf is the probability of finding the atom at a in the state |i a ). Defining an effective current operator, 

l(x 2 ) = 2k" J dx 1 j^(x 1 )D^(x 1 x 2 ), (51) 
we see that the ionization probability reduces to 

i»(fc) = J dx 1 dx 2 SI(x 1 x 2 )(^(x 1 )cb(x 2 )), (52) 
where the spectrometer function x 2 ) is defined in Eq. (|T5l). 



APPENDIX B: MOMENTUM DETECTION PROBABILITY FOR CHARGED PION PAIRS 

In this Appendix we compute the probability for detecting a pair of 7r + : one with momentum k at a and one with 
momentum k' at b. We denote the initial and final states by \T) = \I,i a ,ib,0) and \T) = \k,k' , / a , /bjO), where the 
number of pions in the initial and final states are the same, and k and kl are the measured states. The interaction 
Hamiltonian for the system of detector atom-a and detector atom-b plus two pions is Hj = HJ +H^ + H^, where the 
interactions are given in Appendix A. Two-pion correlation measurements are given by the term fourth-order in Hi 
in the time evolution operator, U(t,t'), which with free fields becomes 

4 f 

- J dx 1 dx 2 dx 3 dx 4 ^(x 1 )T[^( 2 ; 2 )^(x3)]^(x4)T[A AI ( a;i )A,( 2 ; 2 )A CT (x3)A K ( a ;4)]. (53) 

The matrix elements of the pion and photon operators are reduced as follows: We write the pion current-current 
matrix element as 

(k,P\%{x a )%(x 3 )\I) ^idZid^k^'l^ix^ix^Hx^ix^I), (54) 

where d 2 and d% act between the operators carrying the same variables. Using the free-field commutation relation 
[0(x 2 ),0t( X3 )] = V {x 2 X3), a c-number, we express the pion-current operators in terms of the two-pion correlation 
operator plus a term involving T>(x 2 x-s). The matrix element 

(k,k'\^(x 2 )^(x 3 )cf>(x 3 )^x 2 )\I) (55) 
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is exactly the interaction of two pions ionizing two atoms; we remove two pions at two separate space-time points 
and replace them. The commutator term T>(x2 x 3 ) does not contribute to lowest-order scattering, hence we replace 
the matrix element on the right-hand-side of Eq. (|54|) with (|55|); whereupon we eliminate the time-ordering of the 
pion fields in Eq. (m) using the symmetry of these expressions under the interchange X2 «-> £3 and v <-> a. Since we 



require ijj-. at a and Vv, at b, ( |55| ) finally becomes ipUx2)ipl,( x 3)(^\ ( t > ( x 3) c l ) ( x 2)\I}- The vacuum expectation value of 



the time-ordered photon operators in Eq. fl53) factorizes as 
- (OlTiA^A^A^A^Xi)] |0) 



D /JiV (xi,x 2 )D Ka (x4,x 3 ) + D fJ , a (xi,x 3 )D Kl/ {x4:,X2) + D fJiK (xi,x 4: )D l/<7 (x2,x 3 ). (56) 



The last term in Eq. ( p6[) represents photon exchange between detector atoms and is not important for ionization. 

After eliminating the crossed photon lines in the amplitude, and the derivatives of the pion fields as in Appendix A, 
we write 

A^(l — > T) = -e 4 J dx 1 dx 2 dx 3 dx 4 (f Sl \j^(x 1 )\i a }D flv (x 1 x 2 )(/b| 3h \x i )\i h )D KI7 (x i x 3 ) 

x^{x2Wp{x 3 )Ak v k la {Q\<j>{x 3 )4>{x 2 )\I). (57) 

The transition probability is computed in exactly the same way as for single pion detection; redefining the electron 
currents, Eq. (|5l|), and expressing the momentum measurements in terms of spectrometer functions, Eq. (|l5|), we 
derive: 

P^(k,k') = fdx 1 dx2dx 3 dx 4 SI(x 2 x 3 )(^(x 1 ) ( f,\x2)cf > (x 3 )4>(x 4 ))Sl(x 1 x 4 ). (58) 
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